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We study the phase diagram of a two-dimensional ultracold Fermi gas with the synthetic spin- 
orbit coupling (SOC) that has recently been realized at NIST. Due to the coexistence of anisotropic 
SOC and effective Zeeman fields in the NIST scheme, the system shows rich structure of phase sep- 
aration involving exotic gapless superfluid and Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) pairing 
states with different center-of-mass momentum. In particular, we characterize the stability region 
of FFLO states and demonstrate their unique features under SOC. We then show that the effective 
transverse Zeeman field in the NIST scheme can qualitatively change the landscape of the thermo- 
dynamic potential which leads to intriguing effects such as the disappearance of pairing instability, 
the competition between different FFLO states, and the stabilization of a fully gapped FFLO state. 
These interesting features may be probed for example by measuring the in-situ density profiles or 
by the momentum-resolved radio-frequency spectroscopy. 



PACS numbers: 03.75.Ss, 03.75.Lm, 05.30. Fk 

Introduction.- The recent realization of synthetic spin- 
orbit coupling (SOC) in ultracold atoms has greatly ex- 
tended the horizon of quantum simulation in these sys- 
tems [1-3]. Interesting phenomena like topological insu- 
lators, quantum spin Hall effects, topological superfluid- 
ity, etc., where SOC plays a key role [4], may now be 
studied in ultracold Fermi gases. A great amount of ef- 
forts have recently been dedicated to the clarification of 
novel phases and phase transitions in a spin-orbit cou- 
pled, strongly interacting Fermi gas of ultracold atoms, 
where the interplay of SOC, pairing superfluidity and ef- 
fective Zeeman fields can lead to exotic superfluid phases 
in various dimensions [5-22]. 

While most of these studies assume the Rashba SOC, 
only an equal Rashba and Dresselhaus (ERD) SOC has 
been realized in experiments [1-3, 23]. With spatial 
anisotropy and effectively reduced SOC dimension, the 
NIST SOC can lead to qualitatively different superfluid 
phases compared to those under Rashba SOC [21, 22, 24]. 
The picture is further complicated by the existence of 
effective Zeeman fields, both axial and transverse, in 
the NIST scheme, which are tunable by adjusting the 
laser parameters [1]. On the other hand, although there 
has been discussions of Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) state in these systems [19, 20, 24, 25], the FFLO 
state has not been systemtaically characterized in the 
context of the NIST scheme. 

In this work, we study the exotic superfluid phases 
in a two-dimensional (2D) ultracold Fermi gas near a 
wide Feshbach resonance and with the NIST SOC. Given 
the enhanced stability of FFLO states in two dimen- 
sions [26], as well as the recent development of in-situ 
detection schemes [27], 2D Fermi gas is an ideal plat- 
form for the investigation of the unconventional pairing 



states, FFLO in particular, under SOC. By mapping out 
the zero-temperature phase diagram, we find that, in the 
absence of an effective transverse Zeeman field, typically 
two different gapless superfluid states can appear, similar 
to the case of a three-dimensional (3D) polarized Fermi 
gas with Rashba SOC [11, 12]. Furthermore, an FFLO 
state with center-of-mass momentum perpendicular to 
the direction of the anisotropic SOC can be stabilized. 
The quasi-particle (hole) dispersion spectra of this FFLO 
state feature non-trivial gapless contours in momentum 
space, distinct from that of the FFLO state without SOC. 
The unique dispersion spectra and the gapless contours 
may be detected via momentum-resolved radio- frequency 
(rf) spectroscopy [28]. We then discuss the influence of 
effective transverse Zeeman fields on the system. In gen- 
eral, a transverse Zeeman field eradicates the pairing in- 
stability, typical of an attractively interacting Fermi gas 
under SOC in the large polarization limit [17, 29], and 
stabilizes FFLO states with center-of-mass momentum 
along the direction of the anisotropic SOC. These lead 
to the appearance of normal state on the phase diagram, 
and the competition between various FFLO states with 
different center-of-mass momentum. As a result, first- 
order phase transitions are abundant on the phase dia- 
gram, which should leave signatures in the in-situ density 
profiles of a trapped gas. Finally, we identify a continu- 
ous transition between the gapless FFLO states and an 
interesting fully gapped FFLO state, which is the result 
of SOC-induced spin mixing and Fermi surface deforma- 
tion. With recent progress in the experimental investiga- 
tion of 2D Fermi gases [30] and the realization of SOC in 
a degenerate Fermi gas [2] , we expect that many features 
reported in this work will be observed in future experi- 
ments. 
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Model.- We consider a two-dimensional two- 
component Fermi gas with the NIST SOC [1, 2], 
where the Hamiltonian can be written as 

H = ^ $fc«L a ^ + fe(ak t Qkj + h.c.) 
K<?=t,l k 

+ ^2,[(ak x ~ h x )a\^a^ - {ak x - h^a^a^] 

k 

+ 7 E a k+qt a -k+q; a -k'+qi a k'+ q t- (1) 
k,k',q 

Here, £ fc = £fc — A*, £fe = H 2 k 2 /2m, a k(T (a kcr ) is the anni- 
hilation (creation) operator for the hyperfine spin state 
a with a = (t, J,), ni is the atomic mass, S is the quan- 
tization volume in two dimensions, and a denotes the 
strength of the SOC. The effective Zeeman fields h and h x 
are proportional to the effective Rabi- frequency and the 
two-photon detuning, respectively, of the Raman process 
in the experiment [1]. Note that the SOC terms in Hamil- 
tonian (1) differ from the standard ERD form by a spin- 
rotation, and that we adopt the convention in the ERD 
form and label h (h x ) as the axial (transverse) Zeeman 
field. The bare s-wave interaction rate U should be renor- 
malized as [31]: l/U = -l/S£ k l/(.E 6 + 2e k ) ) where E b 
is the binding energy of the two-body bound state in 
two dimensions without SOC, which can be tuned, for 
instance, via the Feshbach resonance technique. 

We focus on the zero-temperature properties of the 
Fulde-Ferrell (FF) pairing states with a single valued 
center-of-mass momentum on the mean-field level [32], 
which should provide a qualitatively correct phase di- 
agram at zero temperature. The effective mean field 
Hamiltonian can then be arranged into a matrix form 

in the hyperfine spin basis |akti a Q 
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where Xj: = ± ak x =p h xi the order parameter Aq = 
^Ek ( a Q-ki a kt)' It i s t nen straightforward to diago- 
nalize the effective Hamiltonian and evaluate the ther- 
modynamic potential at zero temperature 
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where the quasi-particle (hole) dispersion v [y = 1, 2, 
i] = ±) are the eigenvalues of the matrix in Hamiltonian 
(2), and 6(x) is the Heaviside step function. Without 
loss of generality, we assume h,h x > 0, Aq = A, and 






(b) 










V " - 
V 

\ nSFl 




\ 




VAC , 





1 2 

akh/h 



0.5 1 1.5 2 

akh/h 



FIG. 1: (Color online) Phase diagrams on the fi—a plane for 
(a) E b /h = 0.2, h x = 0, and (b) E b /h = 0.5, h x = 0. The solid 
curves are first-order boundaries, while the dashed curves rep- 
resent phase boundaries of continuous phase transitions. The 
dash-dotted curves traversing the nSF phases are the thresh- 
old with A/h — 10 -3 . The axial Zeeman field h is taken to be 
the unit of energy, while the unit of momentum is defined 
through h k%/2m = h. 



Aq to be real throughout the work. The pairing order 
parameter Aq as well as the center-of-mass momentum 
Q for the pairs can then be found by minimizing the 
thermodynamic potential in Eq. (3). 

In general, the Hamiltonian (2) cannot be diagonalized 
analytically, and the thermodynamic potential needs to 
be evaluated numerically. However, for pairing states 
with zero center-of-mass momentum (Q = 0), analytical 
form of the dispersion relation can be obtained along the 
?± 

-E k 2 . Apparently, the branches E^ v can 
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and E 



± 



cross zero for finite A, which leads to gapless superfhiid 
phases with isolated gapless points in momentum space, 
similar to the case of a 3D Fermi gas with Rashba SOC. 
We find numerically that for finite A, the gapless points 
can only occur on the k x = axis [33], therefore it is 
sufficient to study the phase boundaries based on the 
dispersion along this axis. Typically, there can be two 
or four gapless points on the axis, which are symmetric 
with respect to the origin. As the location of the gapless 
points (k yc ) must satisfy the relation EJ^ v = 0, it is 
easy to see that there are four gapless points when /i > 
yjh 2 + h x — A 2 and ti 2 + h\ > A 2 , two gapless points 
v//i 2 + hi - A 2 and h 2 + hi > A 2 . Similar 
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to the case of a 3D Fermi gas with Rashba SOC [12], 
we may then identify nodal superfluid states with two 
gapless points (nSFl) and those with four gapless points 
(nSF2) as distinct gapless superfluid phases. As has been 
pointed out in Ref. [21], these phases have locally non- 
trivial topological properties in momentum space near a 
given gapless point. 

Without effective transverse Zeeman fields- We first 
examine the phases in the absence of effective trans- 
verse Zeeman fields. With h x = 0, the Hamiltonian 
(2) can be diagonalized analytically for pairing states 
with Q = 0. The dispersion spectra are: -E kl = 
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FIG. 2: (Color online) (a) Contour of gapless points in mo- 
mentum space for a standard FFLO state in a polarized Fermi 
gas without SOC. (b) Contours of gapless points in momen- 
tum space for a typical FFLO a state with NIST SOC and with 
h x — 0. (c) Dispersion spectra along the k x — axis for the 
FFLO state in (a), (d) Dispersion spectra along the k x — 
axis for the FFLOy state in (b). The parameters for (b)(d) 
are: E b /h = 0.5, ak h /h = 0.21, fj,/h = 1.83, A Q /h ~ 0.28, 
Qv/kh ~ 0.83, h x = 0. Gapless points on the k x — axis are 
indicated by bold dots. 
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a 2 k 2 x ± 2E , and E^ = -E± 2 , where 
E Q = v/^(/i 2 + « 2 fc2) + A 2 h 2 . While for pairing states 
with finite center-of-mass momentum, one must resort to 
numerical diagonalization. With these, we map out the 
phase diagram on the fi-a plane with fixed h (see Fig. 1). 
Under the local density approximation, the phases tra- 
versed by a downward vertical line in the diagram rep- 
resent those one should encounter starting from a trap 
center and moving to its edge. 

From Fig. 1, we see that the topological superfluid 
phase in a 2D polarized Fermi gas with Rashba SOC is 
now replaced by nSFl and nSF2, as we have anticipated. 
On the other hand, similar to the case of a 2D Fermi gas 
with Rashba SOC, pairing state exists even in the large 
polarization limit [17, 19]. This indicates the persistence 
of pairing instability despite the anisotropy of the NIST 
SOC. Similar to the Rashba case, the pairing instability 
is dictated by the existence of singularities in the gap 
equation when A = [17, 29], which is equivalent to the 
condition that the dispersion spectra should cross zero at 
A = 0, i.e., £| = a 2 k 2 + h 2 . As this is also the condition 
for the existence of Fermi surfaces in a normal state (N) 
under SOC, the normal state is absent from the phase 
diagram and vacuum (VAC) shows up in its place. 

For pairing states with finite center-of-mass momen- 
tum Q, wc find that with h x = 0, Q lies along the y-axis 
(FFLO,,), i.e. perpendicular to the direction of the NIST 
SOC. Interestingly, the pairing states with Q along the 
x-axis (FFLO,) are also stable against the normal state, 




(a) 



gFFLO, 



FFLO, 



N 



VAC 




FIG. 3: (Color online) (a) Phase diagram on the fi—a plane 
for Eb/h = 0.2, h x /h — 0.2. The solid curves are first-order 
boundaries, the dashed curves represent continuous phase 
boundaries. (b)(c) Illustration of dispersion spectra and pos- 
sible pairing states in the presence of h x . While an FFLOa: 
pairing state with a small and negative Q x is possible in (b), 
various FFLO x pairing states with different Q x are possible 
in (c). The thin (solid and dashed) curves represent the non- 
interacting bands formed by hyperfine states under SOC and 
an axial Zeeman field. The bold solid (dashed) curve repre- 
sents the quasi-particle (hole) dispersion. The horizontal line 
represents the Fermi surface. 



but are metastable compared to the FFLO y state. The 
FFLO y state is bounded by the SF state from above and 
the nSF2 state from below, both via first-order phase 
boundaries. While the first-order boundary between the 
FFLO y state and the SF state is consistent with the mean 
field picture in the absence of SOC, the lower boundary 
against the nSF2 state is a direct result of pairing in- 
stability under SOC, as otherwise the nSF2 state at the 
boundary is replaced by the normal state and the phase 
boundary would become continuous. Indeed, the first- 
order boundary between FFLO a and nSF2 approaches a 
critical point at a = 0, which is a second-order critical 
point that can be solved analytically for a 2D polarized 
Fermi gas [26]. As the SOC strength increases, the sta- 
bility region of the FFLO y state decreases on the phase 
diagram, implying competition between FFLOj, pairing 
and SOC. 

To further characterize the properties of the FFLO y 
state, we demonstrate in Fig. 2 typical gapless points 
and dispersion spectra of the FFLO y state in momentum 
space. For a standard FFLO state in a polarized 2D 
Fermi gas, the gapless points form a closed contour on the 
k x -k y plane (see Fig. 2(a)). As SOC breaks the inversion 
symmetry, the dispersion spectra split into four branches 
with the symmetry E^ k v 
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2(b) (d)). Hence two distinct contours of gapless points 
show up in momentum space. This is a unique feature of 
the FFLO y state in a Fermi gas with SOC, which may 
be probed using momentum resolved rf spectroscopy. 

With an effective transverse Zeeman field- We now 
turn to the case where an effective transverse Zeeman 
field h x is present. As illustrated in Fig. 3(b) (c), the 
transverse field makes the bands asymmetric. The Fermi 
surface becomes deformed, which may render s-wave 
pairing with zero center-of-mass momentum unfavorable 
even under SOC, i.e., pairing instability may no longer 
exist. This is reflected in the thermodynamic poten- 
tial, where a local minimum can appear at A = 0. A 
more thorough examination shows an intriguing land- 
scape, where the thermodynamic potential may have up 
to three local minima with Q = 0, corresponding to the 
SF, the nSF2 and the normal state, respectively [34]. 
Hence, there can potentially be three different first-order 
transitions with Q = 0. 

The picture is further complicated when we take FF 
pairing states into account. With deformed Fermi sur- 
face, pairing states with center-of-mass momentum along 
the x-direction can be stabilized against the SF state. A 
qualitative understanding is that due to the Fermi sur- 
face asymmetry, the local minima in the case of Q = 
can be shifted to finite Q x in the presence of a finite h x . 
The magnitude and direction of this shift depends on 
the Fermi surface asymmetry (see Fig. 3(b)(c)). Under 
certain conditions, there may even exist several different 
FFLOz states with distinct Q x and Aq (see Fig. 3(c)). 
The ground state of the system therefore is the result of 
competitions between the various FFLO pairing states 
and the normal state. 

With this understanding, we calculate the typical 
phase diagram for finite h x (see Fig. 3(a)). Apparently, 
the most striking effect of the transverse field is the sta- 
bilization of FFLOa, states over a large parameter region, 
with the center-of-mass momentum of the FFLO^ states 
opposite to the direction of the N1ST SOC (Q x < 0). 
This is qualitatively consistent with the findings in Ref. 
[20], where a 3D Fermi gas with Rashba SOC and trans- 
verse fields was considered. When the SOC strength in- 
creases, the magnitude of Q x decreases exponentially as 
the Fermi surface asymmetry decreases, and the energy 
of the FFLOa; state becomes exponentially close to that 
of an SF state (see Fig. 4(a)(b)). Similar scaling re- 
lations can be found with increasing \i [34]. As a or // 
increases, when the magnitude of Q x becomes sufficiently 
small, the FFLO x state undergoes a continuous transi- 
tion and becomes fully gapped. In the weak coupling 
limit, this gapped FFLO state (gFFLO^) can be under- 
stood as a pairing state within one helicity band, where 
the Fermi surface deformation induced by transverse field 
is accommodated by a finite center-of-mass momentum 
(see Fig. 3(b)). We expect that such an exotic pair- 
ing state is characteristic of the co-existence of SOC and 
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FIG. 4: (Color online) (a) Evolution of the center-of-mass mo- 
mentum and pairing order parameter (inset) for the ground 
state as the SOC strength increases, (b) Evolution of the 
difference in thermodynamic potential between the SF state 
and the FFLOa; state with increasing SOC strength, with 
SQ = SIsf — UfflOx- The inset shows the evolution of 
the minimum excitation gap with increasing SOC strength. 
The parameters for (a)(b) are: Et/h — 0.2, h x /h = 0.2, 
(i = 0. (c) Typical gapless contours in momentum space for 
an FFLO B state in the presence of transverse Zeeman field, 
with the parameters: Eb/h = 0.5, h x /h — 0.5, akh/h — 0.4, 
n/h = 2.22, A Q /h ~ 0.32, Q v /k h ~ 0.86. (d) Typical gap- 
less contours in momentum space for an FFLOa; state in the 
presence of transverse Zeeman field, with the parameters: 
E b /h = 0.2, h x /h = 0.2, ak h /h = 0.9, fi/h = 2, A Q /h ~ 0.19, 
Q x /k h ~ -0.11. 



Fermi surface asymmetry, and is not limited to systems 
under NIST SOC. On the other hand, the existence of 
various first-order boundaries and end points in Fig. 3(a) 
suggests stabilization of FFLO states with different Q x 
and Aq, consistent with our previous analysis. These 
FFLOz states eventually merge into a single FFLO^ state 
beyond the end points. For finite h x , we can still find 
considerable stability region for the FFLOj, state, which 
is bounded from the FFLO^ or gFFLO^ states by first- 
order boundaries. The typical dispersion of the FFLO y 
state here has the symmetry E^ jkyjV = -E v l k ^ Qv _ ky y 
This is reflected in the contours of the gapless points (see 
Fig. 4(c)), which is qualitatively different from that of 
an FFLO x state (see Fig. 4(d)). 

Conclusion- We have studied the exotic phases in a 
2D Fermi gas under the NIST SOC. The competition 
between the effective Zeeman fields and the anisotropic 
SOC leads to complicated phase structures, implying 
multiple phase separated states for a uniform gas and 
rich spatial phase structure for a trapped gas. Impor- 
tantly, we demonstrate that various FFLO states, includ- 
ing a fully gapped FFLO state, can be stabilized under 
the NIST SOC. While the unique dispersion spectra and 
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gapless contours of the FFLO states may serve as signa- 
tures for experimental detection, the various first-order 
phase transitions should also leave signatures in the in- 
situ density profiles of a trapped gas. 

Acknowledgements.- We thank Lin Dong, Han Pu, and 
Vijay B. Shenoy for helpful discussions. This work is 
supported by NFRP (2011CB921200, 2011CBA00200), 
NNSF (60921091), NKBRP (2013CB922000), NSFC 
(11105134, 11274009), the Fundamental Research 
Funds for the Central Universities (WK2470000001, 
WK2470000006), and the Research Funds of Renmin 
University of China (10XNL016). W.Z. would also like 
to thank the NCET program for support. 



* Electronic address: wzhangl@ruc.edu.cn 
' Electronic address: wyiz@ustc.edu.cn 
[1] Y.-J. Lin, K. Jimenez-Garcfa, and I. B. Spielman, Nature 

(London) 471, 83 (2011). 
[2] P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. 
Zhai, and J. Zhang, Phys. Rev. Lett. 109, 095301 (2012); 
L. W. Cheuk, A. T. Sommer, Z. Hadzibabic, T. Yef- 
sah, W. S. Bakr, and M. W. Zwierlein, ibid 109, 095302 
(2012). 

[3] J.-Y. Zhang, et al, Phys. Rev. Lett. 109, 115301 (2012). 
[4] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010) . 

[5] C. Zhang, S. Tewari, R. M. Lutchyn, and S. Das Sarma, 

Phys. Rev. Lett. 101, 160401 (2008). 
[6] M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 

103, 020401 (2009). 
[7] J. P. Vyasanakere, S. Zhang, and V. B. Shenoy, Phys. 

Rev. B 84, 014512 (2011). 
[8] M. Gong, S. Tewari, and C. Zhang, Phys. Rev. Lett. 107, 

195303 (2011). 

[9] Z.-Q. Yu and H. Zhai, Phys. Rev. Lett. 107, 195305 

(2011) . 

[10] H. Hu, L. Jiang, X.-J. Liu, and H. Pu, Phys. Rev. Lett. 

107, 195304 (2011). 
[11] M. Iskin and A. L. Subasi, Phys. Rev. Lett. 107, 050402 

(2011). 

[12] W. Yi and G.-C. Guo, Phys. Rev. A 84, 031608(R) 

(2011) . 

[13] L. DellAnna, G. Mazzarella, and L. Salasnich, Phys. 

Rev. A 84, 033633 (2011). 
[14] M. Gong, G. Chen, S. Jia, and C. Zhang, Phys. Rev. 

Lett. 109, 105302 (2012). 
[15] L. Han and C. A. R. Sa de Melo, Phys. Rev. A 85 

011606(R) (2012). 
[16] L. He and X.-G. Huang, Phys. Rev. Lett. 108, 145302 

(2012) . 

[17] J. Zhou, W. Zhang, W. Yi, Phys. Rev. A 84, 063603 

(2011) . 

[18] X. Yang and S. Wan, Phys. Rev. A 85, 023633 (2012). 
[19] W. Yi and W. Zhang, Phys. Rev. Lett. 109, 140402 

(2012) . 

[20] Z. Zheng, M. Gong, X. Zou, C. Zhang, and G.-C. Guo, 

arXiv:1208.2029. 
[21] L. Han and C. A. R. Sa de Melo, arXiv: 1206.4984. 
[22] M. Iskin and A. L. Subasi, arXiv: 121 1.4020. 



[23] X.-J. Liu, M. F. Borunda, X. Liu, and J. Sinova, Phys. 

Rev. Lett. 102, 046402 (2009). 
[24] V. B. Shenoy, arXiv:1211.1831. 

[25] L. Dong, L. Jiang, H. Hu, and H. Pu, arXiv:1211.1700. 
[26] G. J. Conduit, P. H. Conlon, and B. D. Simons, Phys. 

Rev. A 77, 053617 (2008). 
[27] N. Gemelke, X. Zhang, C.-L. Hung, and C. Chin, Nature 

460, 995 (2009). 
[28] J. T. Stewart, J. P. Gaebler, and D. S. Jin, Nature 454, 

744 (2008). 

[29] S. Tewari, T. D. Stanescu, J. D. Sau, and S. Das Sarma, 

New J. Phys. 13, 065004 (2011). 
[30] K. Martiyanov, V. Makhalov, and A. Turlapov, Phys. 

Rev. Lett. 105, 030404 (2010); B. Frohlich, M. Feld, E. 

Vogt, M. Koschorreck, W. Zwerger, and M. Kohl, ibid 

106, 105301 (2011); P. Dyke, E. D. Kuhnle, S. Whitlock, 

H. Hu, M. Mark, S. Hoinka, M. Lingham, P. Hannaford, 
and C. J. Vale, ibid 106, 105304 (2011); A. T. Som- 
mer, L. W. Cheuk, M. J.-H. Ku, W. S. Bakr, and M. W. 
Zwierlein, ibid 108, 045302 (2012); Y. Zhang, W. Ong, 

I. Arakelyan, J. E. Thomas, ibid 108, 235302 (2012). 
[31] M. Randeria, J.-M. Duan, and L.-Y. Shieh, Phys. Rev. 

Lett. 62, 981 (1989). 
[32] P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964); 
A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 
20, 762 (1965). 

[33] For the case of zero center-of-mass pairing and with no 
transverse field, it can be shown directly from the an- 
alytic form of the dispersion spectra that the gapless 
points can only occur on the k x = axis. 

[34] See Supplemental Material for details. 



G 




FIG. 5: (Color online) (a) Phase diagram on the fi-a plane for Eb/h = 0.2, hx/h — 0.2. Only zero center-of-mass momen- 
tum pairing states are considered here. The solid curves are first-order boundaries, while the dashed curves represent phase 
boundaries of continuous phase transitions, (b) Typical thermodynamic potentials with up to three local minima, with the 
parameters: akh/h — 0.33, fi/h = 4.53 (solid curve); akh/h = 0.32, \i/h — 4.53 (dotted curve); akh/h = 0.42, fi/h = 4.2 
(dash-dotted curve); akh/h = 0.33, \i/h = 4.6 (dashed curve); while Eb/h = 0.2, h x /h = 0.2 for all cases, (c) Thermodynamic 
potential for Q — and Eb/h = 0.2, h x /h — 0.2, akh/h ~ 2.12, (i/h ~ —0.98. The local minima are located at A/h ~ 0.3 and 
A/h = 0. (d) Thermodynamic potential as a function of Aq and Q x for Eb/h = 0.2, h x /h = 0.2, akh/h ~ 2.12, jt*//i ~ —0.98. 
The local minima are located at (Aq/A ~ 0.29, Q x /kh ~ —0.11) and (Aq//j ~ 0.16, Q x /kh ~ —0.35). 



Supplementary material 

In this supplementary material, we provide more details regarding the qualitative understanding and the properties 
of the FFLOz state under the NIST scheme. 

As we have discussed in the main text, the various FFLO^ states can be qualitatively understood as shifts of the 
local minima to finite Q x induced by Fermi surface asymmetries along the x-axis. This can be more easily understood 
by comparing the Q = phase diagram (see Fig. 5(a)) with the phase diagram in the main text (Fig. 3(a) therein). 
In the presence of a transverse Zeeman field, up to three local minima can be found in the thermodynamic potential 
for Q — 0, corresponding to the normal, the nSF2 and the SF states, respectively (see Fig. 5(b)). This is reflected 
in the various first-order boundaries in Fig. 5(a). When we consider finite center-of-mass momentum pairing, these 
local minima will be shifted into the plane of finite Q x . As a result, the first-order boundaries in Fig. 3(a) of the 
main text are reminiscent of the first-order boundaries for Q = in Fig. 5(a), but are shifted slightly. Of course, the 
phase diagram would be further modified by the competition from the FFLOj, states, and the various local minima 
along Q = Q x x may merge as the parameters are tuned, which generates the end points we see in Fig. 3(a) of the 
main text. As an example, we plot in Fig. 5(c) a typical thermodynamic potential with Q = and in Fig. 5(d) the 
corresponding thermodynamic potential landscape when finite Q x pairing is considered. It appears that the original 
two local minima are both shifted into the finite Q x plane. 

A remarkable feature of the phase diagram in the presence of a transverse Zeeman field and SOC is the enhanced 
stability region of the FFLOz state which seems to replace the SF state in the Q = phase diagram. As we have 
discussed in the main text, this is due to the Fermi surface deformation. As the Fermi surface deformation becomes 
smaller with increasing chemical potential or with increasing SOC strength, it is natural to expect that the energy 
gain by adopting a finite center-of-mass momentum pairing in the FFLO^ state should become less. This can be 
seen by a direct scaling of the pairing parameters Q x and Aq of the FFLO^ state with increasing SOC strength 
(see Fig. 4(a) in the main text) or with increasing chemical potential (see Fig. 6(a)). Apparently, while the pairing 
order parameter Aq keeps increasing with chemical potential and SOC strength, the magnitude of the center-of-mass 
momentum of the FFLOa; state becomes exponentially small in the large fj, and/or large a limit. Furthermore, the 
difference in thermodynamic potential between the SF state and the FFLOz state also becomes exponentially small. 
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FIG. 6: (Color online) (a) Evolution of the center-of-mass momentum and the pairing order parameter (inset) for the ground 
state as the chemical potential increases, (b) Evolution of the difference in thermodynamic potential between the SF state and 
the FFLO:r state with increasing chemical potential, with 8Q = S7sf — ^fflOx- The parameters are: Eb/h = 0.2, h x /h = 0.2, 
and akh/h = 1.9. 

Therefore in real experiments, the parameter region where the FFLO^ states can be observed should feature moderate 
SOC and moderate chemical potential. 



